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ABSTRACT 

Heat  and  mass  transfer  characteristics  of  water  based 
nanofluids  over  a  nonlinear  shrinking  surface  with 
suction  is  studied  numerically  in  the  incidence  of 
magnetic  field,  thermal  radiation  under  the  boundary 
conditions  heat  flux  and  mass  flux.  Similarity 
transformations  are  used  to  renovate  the  governing 
partial  differential  equations  into  higher  order 
nonlinear  ordinary  differential  equations.  These 
equations  are  numerically  solved  using  fourth  order 
RungeKutta  technique  along  with  Nachistem  - 
Swigert  shooting  method.  Numerical  results  are 
achieved  for  the  nanofluid  velocity,  nanofluid 
temperature  and  nanoparticle  concentration  as  well  as 
for  the  skin  friction  coefficient,  wall  temperature  and 
concentration  of  the  wall  for  various  values  of 
governing  physical  parameters  viz.  magnetic 
interaction  parameter,  suction  parameter,  shrinking 
parameter,  thermal  radiation  parameter,  temperature 
ratio  parameter  and  solid  volume  fraction  of  the 
nanoparticles  such  as  copper  and  silver. 

Keywords:  Nanofluid,  MHD,  Radiation,  Shrinking 
sheet,  heat  flux,  mass  flux 

I.  INTRODUCTION 

Research  of  boundary  layer  flow  over  a  shrinking 
sheet  has  very  much  interest  on  last  few  years  due  to 
its  vast  range  of  practical  applications  in  industries. 
Shrinking  sheet  is  a  surface  which  decreases  in  size  to 
a  certain  area  due  to  an  imposed  suction  or  external 
heat.  One  of  the  most  common  applications  of 
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shrinking  sheet  problems  in  industries  is  shrinking 
flim.  Shrinking  flim  phenomenon  is  very  supportive 
in  packaging  bulk  products  as  it  can  be  unwrapped 
easily  with  adequate  heat.  The  viscous  boundary  layer 
flow  over  a  shrinking  sheet  was  first  coined  by  wang 
(1990)[1].  In  this  work,  he  concluded  that  the 
boundary  layer  flow  over  a  shrinking  surface  have 
suction  effect  on  the  boundary.  Following  this 
pioneering  work,  many  researchers  have  been 
investigated  on  the  research  area  of  shrinking  sheet. 
Some  of  the  researchers  are  Miklav  et.al  (2006)  [2], 
wang  (2008)  [3],  Fang  (2009)  [4],  Arifin  et.al  (2011) 
[5],  Rohniet.al(2012)  [6],  Bachok  et.  al  (2012)  [7], 
Anuarlshak  et.  al  (2013)  [8],  Kalidas  Das  (2014)  [9] 
and  Rahman  et  al.(2016)  [10], 

Magneto  hydrodynamics  (MHD)  has  many 
engineering  applications  in  the  real  world  such  as 
astrophysics(planetary  magnetic  field),  MHD  pumps, 
MHD  generators,  MHD  flow  meters,  jet  printers, 
fusion  reactors  and  so  on.  As  a  result  of  the  wide 
applications  of  MHD,  MHD  boundary  layer  flow  of 
viscous  fluid  over  a  stretching  sheet  was  studied  by 
numerous  researchers.  Andersson  et  al.(1992)  [11] 
studied  the  effect  of  hydromagnetic  flow  of  a  power 
law  fluid  over  a  stretching  sheet.  Hydromagnetic 
boundary  layer  flow  of  a  viscous  incompressible  fluid 
which  is  caused  by  a  sheet  stretching  according  to  a 
power  law  velocity  distribution  in  the  presence  of  a 
magnetic  field  was  investigated  by  Chaim  (1995) 
[12],  Thermal  radiation  effects  on  an  electrically 
conduction  fluid  arise  in  many  physical  or  practical 
applications  such  as  electrical  power  generation,  solar 


@  IJTSRD  |  Available  Online  @  www.ijtsrd.com  |  Volume  -  1  |  Issue  -  6  |  Sep  -  Oct  2017 


Page:  1302 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 


power  technology,  nuclear  reactors  and  nuclear  water 
disposed.  Radiation  effects  on  boundary  layer  flow 
with  or  without  considering  the  magnetic  field  have 
been  researched  by  several  authors  under  certain 
contions.  The  hydromagnetic  flow  fast  a  flat  plate  by 
the  presence  of  radiation  was  analysed  by  Raptis  and 
Massallas  (1998)  [13], The  problem  of  radiation 
effects  on  the  heat  transfer  over  a  stretching  surface 
has  been  solved  by  Elbashbeshy  (2000)  [14].  Heat 
transfer  analysis  for  fluid  flow  over  an  exponentially 
stretching  porous  sheet  with  surface  heat  flux  in 
porous  medium  was  worked  by  Iswar  Chandra 
Mandal  (2013)  [15].  Maiti  (2014)  [16]  was  published 
the  effect  of  radiative  heat  flux  of  a  magneto 
hydrodynamics  micropolar  fluid  flow  towards  a 
stretching  or  shrinking  vertical  surface  in  presence  of 
a  heat  source  or  sink.  Stagination  point  flow  of  an 
electrically  conducting  micropolar  fluid  with  MHD 
and  radiation  towads  a  horizontal  shrinking  sheet  was 
published  by  Das  and  Sarkar  (2016)  [17]. 

A  word  nano  fluid  was  first  defined  by  Choi  (1995) 
[18].  Nano  fluid  is  a  new  domain  of  heat  transfer  fluids 
containing  nanometer  sized  (1  -  lOOnm)  particles  that 
are  stable  suspended  in  the  conventional  fluids  like 
water,  oil,  ethylene  glycol.  Commonly  used 
nanoparticles  are  oxide  ceramics  (Al203,Cu0), metals 
(Cu,Ag,Au),  Semi  conductors  ((Ti02,SiC),  nitride 
ceramics(AlN,SiN)  and  etc.  Due  to  nanofluids  novel 
properties,  it  can  be  widely  used  for  various  heat 
transfer  applications  such  as  solar  energy  will 
maximized  with  the  change  of  nanoparticle 
size, shape,  volume  fraction  of  nanoparticles, 
automotive  and  air  conditioning  cooling,  cooling  of 
transfer  oil,  improving  diesel  generator  efficieny  in 
nuclear  reactor,  defense  and  space.  The  research  about 
nanofluid  was  reported  by  several  researchers 
rYiminXuan  (2000)  [19],  Das  et.  al(2003)  [20], Phelan 
et.  al  (2005)  [21], Ding  et.  al  (2007)  [22], 
Manca(2010)  [23]].  The  unsteady  stagination  point 
flow  over  a  shrinking/stretching  sheet  with  heat 
transfer  in  the  presence  of  prescribed  heat  flux  was 
researched  by  Sualiet.  al  (2012)  [24],  Winifred 
NdukuMutuku  (2014)  [25]  published  the  effects  of 
MHD  nonlinear  boundary  layer  flow  and  heat  transfer 
of  nanofluids  past  a  permeable  moving  flat  surface 
with  thermal  radiation  and  viscous  dissipation.  The 
problem  of  MHD  boundary  layer  flow  of  copper 
water  nanofluid  and  alumina  water  nanofluid  over  a 
moving  permeable  flat  plate  has  been  solved  by 
Eshetuhaile  and  Shankar  (2015)  [26], 

SandeepNaramgar  and  Sulochana  (2016)  [27] 


investigated  about  the  nanofluid  flow  past  a  stretching 
and  shrinking  sheet  under  the  consideration  of 
magnetic  field  with  suction/injection. 

Hence  there  is  no  work  has  been  researched  on  steady, 
two  dimensional  flow  of  nanofluids  over  a  nonlinear 
shrinking  sheet  with  nonlinear  thermal  radiation  and 
MHD  in  the  presence  of  heat  flux  and  mass  flux.  This 
study  focused  on  nanofluids  is  the  water  based  metal 
nanoparticles  such  as  Cu  and  Ag. 

II.  MATHEMATICAL  FORMULATIONS 
A.  Governing  Equations 

The  problem  of  steady,  hydromagnetic 
incompressible,  viscous,  two  dimensional  boundary 
layer  flow  of  a  nanofluid  over  a  nonlinear  shrinking 
sheet  with  prescribed  heat  flux  and  mass  flux  is 
considered.  The  x-  axis  chosen  in  the  direction  of  the 
flow  along  the  nonlinear  shrinking  surface  and  y-  axis 
is  chosen  normal  to  the  shrinking  surface  with 
velocity  components  u  =  uw(x)  =  cxn,  c  >  0  and  v  in 
these  directions. The  nonlinear  shrinking  surface  is 
subjected  to  suction.  Two  types  of  water  based 
nanofluids  are  considered  for  this  work.  The 
Schematic  diagram  of  the  problem  is  given  in  Fig.  (1). 


Assume  that  the  non  uniform  transverse  magnetic 

n- 1 

field  of  strength B  ^  =  B°x  2  (Chaim  (1995))),  B0  is 
the  uniform  magnetic  field  is  applied  in  the  y  - 
direction.  It  is  also  assumed  that  the  ambient 

T 

temperature  is  constant,  denoted  by  00  and 

rji  _  rp  1  m 

»  "  »  is  the  surface  temperature  with  m  being 
surface  temperature  parameter  and  b  is  a  positive 
constant.  Consider  the  nanofluids  to  be  gray, 
absorbing  and  emitting  radiation  but  non  scattering 
medium.  Also  the  radiative  heat  flux  qr  is  negligible  in 
x  axis  when  compared  to  that  in  the  y  direction  in  the 
energy  equation  is  described  using  Rosseland 
approximation. 
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Hence,  under  the  above  consideration  and  boundary 
layer  assumptions,  the  governing  equations  of  the 
hydromagnetic  flow  of  nanofluid  over  a  nonlinear 
shrinking  surface  is  given  by 


Following  the  lines  of  Hady  et  al.  (2012)  [28]  the 
effective  density,  specific  heat  capacity,  dynamic 
viscosity  &  thermal  conductivity  of  the  nanofluid  are 
given  by 
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where  u 

and 

v  are  the 

velocity  components 

and  y  axes  respectively,  a  is  the  electrical 

conductivity,  '“vis  the  effective  dynamic  viscosity  of 
the  nanofluid,  pnf  is  the  density  of  the  nanofluid, anf  is 
the  thermal  diffusivity  of  the  nanofluid,  cp  is  the 
specific  heat  at  constant  pressure,  D  is  the 
thermodiffusion  coefficient.  Using  Rosseland 
approximation,  the  net  radiative  heat  flux  is  written 
as, 


16  <r  *  (  T  3  8JT  j 

3  k*  l  dy  J 


16  a 


3  k* 


3  T 


d2  T 


(5) 


Where  a*  is  the  Stefan-Boltzmann  constant  and  k*  is 
the  mean  absorption  coefficient. 


B.  Thermophysical  properties  of  the  Nanofluid 

The  present  study,  water  is  considered  as  the  base 
fluid  with  Prandtl  number  Pr  =  6.2.  The  nanofluids  are 
the  mixer  of  water  with  solid  spherical  copper  and 
silver  nanoparticles. 

The  viscosity,  specific  heat  capacity,  density,  thermal 
diffusivity  and  thermal  conductivity  of  the  nanofluids 
depend  on  the  volume  fraction  cp  of  the  nanoparticles 
used. 


The  thermal  diffusivity  of  the  nanofluid  is  given  by 

k«  (6) 


-  'if 


(Pc„)„f 

Where  Ay/and  (pcp)nf  are  the  thermal  conductivity  of 
the  nanofluid  and  specific  heat  capacity  of  the 
nanofluid 


Pnf  =  0  -  <t>)Pf  +  ‘PPs  ’ 

(Pcp  )„/  =  (1  -  <f>)(PCp  )  f  +  <Kpcp  )s 

ft,  k  =k  (*,  +  2*/)-2^(fc/-fr,) 

Mnf  =  (1  -  <02S  ’  nf  f  (ks  +  2kf)  +  </>(kf-ks) 

(V) 

Where,  p  is  the  density  of  the  nanofluid,  p  is  the 
effective  viscosity  of  the  nanofluid,  vnf  is  the 

kinematic  viscosity  of  the  nanofluid,  (p  is  the  solid 
volume  fraction,  cp  is  the  specific  heat  at  constant 
pressure,  p  and  p  are  the  densities  of  the  base  fluid 

and  the  nanoparticle.  ( pcp)f  and  {Pcp)s  are  the 

specific  heat  capacities  of  the  base  fluid  and  the 
nanoparticles  respectively  and  kf  and  ks  are  the 
thermal  conductivities  of  the  base  fluid  and  the 
nanoparticles  respectively.  The  relation  for  effective 
dynamic  viscosity  Pnf  used  in  the  present  study  is 

taken  from  the  equation  by  Brinkman  (1952).  The 
Thermophysical  properties  of  water  and  nanoparticles 
are  tabulated  in  Table  1 . 


Table  1:  Thermo-Physical  Properties  Of  Fluid  And 
Nanoparticlesat  25°C 


Physical 

properties 

Water 

fluid 

Cu 

Ag 

cP 

4179 

385 

235 

P 

997.1 

8933 

10500 

K 

0.613 

400 

429 

a 

1.47 

1163.1 

1738.6 

C.  Boundary  conditions 

The  boundary  conditions  of  the  prescribed  heat  flux 
and  mass  flux  of  the  problem  is  defined  as  follows, 

-k  —  =  a  And  _  n  -  m  ,  here  surface 

nf  dy  y"  °  dy  ~  w 

heat  flux  qw  and  surface  mass  flux  mw. 

Therefore,  the  appropriate  boundary  conditions  are, 

u  =  ~uw(x)  ’V  =  -r„(x)  ’  -knf  ^3L=  q„  ’  - D  =  mw 

dy  dy 

at  y  =  0  (8) 

a  — y  0 ,  T  — y  1,  ,  C  ^  C  -  /  a  s  y  — >  oo 
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M*) 

VA*)  = 


is  the  variable 


C(n  +  \)v  ,  izL 


suction  velocity. 


D.  Similarity  transformations 

To  assist  the  analysis,  the  following  similarity 
transformations  are  introduced. 
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077T); 


fin) 
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— — — x  2  dir/) 
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2v  — 

f  x  2  (f>ir/)  (9) 


The  horizontal  and  vertical  velocity  components  are 
derived  by  using  stream  functions 


u  =  ffff—  =  Cxn  f  (r/)  and 


di/r  _  fc(n  +  \)v 


2-x  2 


f(v)  +  - — 

n  +  1 


(10) 


The  continuity  equation  is  satisfied  identically. 
Therefore,  the  transformed  momentum  equation, 
energy  equation  and  concentration  equation  is  known 
by 

v  '  (11) 


R  rtf  J  R  nf 


+  Yr-f(<t,ffe'-2pfe))  =  0 


i"  +  Scff  =  0 
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(13) 


Where 


=  i  -  <t>  +  <t> 
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with  the  corresponding  boundary  conditions  obtain 
from  equation  (5)  in  the  following  form 


/(0)=  /„>/'(O)=-l>0'(O)~-W'(O)=  -1  at 
r/  =0, 

/  '  (oo  )  =  o ,  ^  (co  )  — >  o  (^(oo)— >  0  as  r/  ^  go 


The  non  dimensional  parameters  appeared  in 
equations  (11)  to  (14)  are  defined  as  follows 

2  _  2aBl  js  the  magnetic  interaction 
c(n  +  l)p, 

parameter,  N  =  k/k  *  is  the  radiation  parameter 

*  4a  *  rj 

Pr  =  v  /  ipc p )/  js  the  Prandtl  number  _  Tw  is  the 

k  uw 

temperature  ratio  parameter,  ^  =  2  n  is  the 

shrinking  parameter, fw  is  the  suction  parameter, 
„  i/  ,  is  the  Schmidt  number. 

Sc  —  ~ 

D 

E.  Physical  quantities  of  interest 

The  quantities  of  engineering  interest  in  this  study  are 
the  non-dimensional  skin  friction  coefficient 
r  _  where  ( 8u)  (15) 

f~  Pf<  T^^[e7L 

C/K)1/2  =  (S/"(°)’  the 

Reynolds  number 

III.  NUMERICAL  SOLUTIONS 

The  set  of  nonlinear  coupled  ordinary  differential 
equations  (11)  to  (13)  subject  to  the  boundary 
conditions  in  (14)  constitute  a  nonlinear  boundary 
value  problem  which  is  reduced  to  an  initial  value 
problem  to  start  with.  In  order  to  solve  these 
numerically,  the  most  efficient  numerical  shooting 
technique  such  as  the  Nachtshiem-Swigert  iteration 
shooting  scheme  is  utilized  to  satisfy  the  asymptotic 
boundary  conditions  at  r/  — >co .  The  appropriate  finite 
values  of  r/  — »oo  are  chosen.  Some  initial  conditions 
are  required  to  solve  the  above  equations.  Hence  three 
initial  conditions  /"( 0)  ,  0(0)  ,  (f)(0)  which  are  not 
prescribed  are  to  be  obtained  using  Nachtshiem- 
Swigert  shooting  iteration  technique.  Now  the 
boundary  value  problem  is  converted  into  initial  value 
problem  and  has  been  solved  numerically  using  fourth 
order  Runge-Kutta  integration  method.  The  crux  of 
the  problem  lies  in  making  good  guesses  for  /"( 0)  , 
0(0)  ,  (p(0) .  They  are  guessed  properly  to  get  the 
convergence  and  solution  for  various  values  of 
physical  parameters.  The  accuracy  level  is  chosen  up 
to  1 0~5 .  | 
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IV.  RESULTS  AND  DISCUSSIONS 

The  accuracy  of  the  present  method  was  authorized  by 
comparisons  with  numerical  values  of  local  Nusselt 
number  by  previously  published  articles  Ishak  et 
al.(2008)  [29]and  Elsayed  et  al.(2010)[30]  in  the 
absence  of  volume  fraction,  without  taking  of 
radiation  and  the  value  of  temperature  ratio  parameter 
is  1.0  and  the  shrinking  parameter  is  14.  The  numerical 
results  are  in  Table  2.  From  this  table,  our  results  are 
good  concurrence  with  the  authors’  results  Ishak  et  al. 
and  Elsayed  et  al.  for  different  values  of  Prandtl 


Pr 

Ishak  et 
al.  [29] 

Elsayed  et 
al.[30] 

Present 

results 

0.72 

0.8086 

0.808 

0.80863 

1.0 

1.0 

1.0 

1.0 

10.0 

3.7202 

3.7207 

3.72067 

number 


TABLE  2:  COMPARISON  OF  LOCAL 
NUSSELT  NUMBER  AT  c()=0.0,  NR->oo,  0W  =  1.0 
AND  p  =  i  AT  DIFFERENT  VALUES  OF  PR 
WITHPRE VIOU SLY  RESEARCHED  DATA 

The  numerical  solution  of  the  problem  is  found  and 
the  effects  of  the  physical  parameters  on  the  flow, 
temperature  and  concentration  are  depicted  through 
graphs.  Also  the  engineering  interest  of  physical 
quantities  such  as  skin  friction  coefficient,  wall 
temperature  and  concentration  of  the  wall  is  given  in 
tabular  form. 

In  this  problem,  the  nanoparticles  solid  volume 
fraction  is  fixed  in  the  range  of  0<(j)<0.2,  magnetic 
interaction  parameter  M2  =  0.5,  1.0, 2.0,  2. 5, shrinking 
parameter  p  =  -1.0,  1.0,  1.5, 2.0,  suction  parameter  fw 
=  2.0,  2.3,  2.4,  2.5,  The  radiation  parameter  Nr  =  1.0, 
1.5,  2.0,  2.5,  The  temperature  ratio  parameter  0W  = 
0.8,  0.9,  1.1,  1.2  and  Schmidt  number  Sc  =  0.3,  0.5, 
0.7,  0.8  are  obtained  and  demonstrated  by  means  of 
graphs. 

The  effect  of  magnetic  field  on  velocity,  temperature 
and  concentration  of  nanofluids  for  both  copper  water 
nanofluid  and  silver  water  nanofluid  is  shown  in 
figs.2,  3  and  4.  It  is  noted  that  in  the  shrinking  surface 
the  nanofluid  velocity  increased  for  increasing 
magnetic  parameter.  From  fig.3  and  fig.4  the 
nanofluids  temperature  and  concentration  decreases 
due  to  the  effect  of  magnetic  field.  But  the  thermal 


boundary  layer  thickness  has  very  small  amount  of 
variation  while  increasing  the  values  of  magnetic 
interaction  parameter. 

The  impact  of  shrinking  parameter  over  dimensionless 
velocity  is  portrayed  in  fig.5.  The  figure  proves  that 
the  nonlinear  shrinking  parameter  shrinks  the  velocity 
of  both  copper  water  nanofluid  and  silver  water 
nanofluid.  Fig.6  illustrates  the  dimensionless 
temperature  distribution  for  various  values  of 
shrinking  parameter.  For  increased  values  of 
shrinking  parameter,  the  nanofluid  temperature  lifts 
and  the  thermal  boundary  layer  thickness  increased. 
The  effects  of  shrinking  parameter  on  nanoparticle 
concentration  are  displayed  in  fig.7.  For  both 
nanofluids,  the  shrinking  parameter  rises  the 
nanoparticle  concentration  when  <))  =  0.05. 

Fig.8  highlights  the  effects  of  suction  parameter  on 
the  dimensionless  velocity  distribution.  In  the 
shrinking  surface  the  values  of  suction  parameter 
augments  the  nanofluid  velocity  for  both  nanofluids. 
Variation  in  dimensionless  temperature  and 
concentration  distribution  is  visualized  through  figs.  9 
and  10.  For  enhances  the  values  of  suction  parameter 
(fw>  0),  the  thermal  boundary  layer  thickness  and 
concentration  of  both  copper  water  nanofluid  and 
silver  water  nanofluid  gets  decelerated. 

Figs.  11,  12  and  13  elucidates  the  impact  of  solid 
volume  fraction  over  the  dimensionless  velocity, 
temperature  and  concentration  of  copper  water 
nanofluid  and  silver  water  nanofluid.  The  nanofluid 
velocity  enlarges  due  to  the  variation  in  solid  volume 
fraction  is  shown  in  fig.ll.  The  thermal  boundary 
layer  thickness  expands  when  the  solid  volume 
fraction  enhance.  In  fig.  13  The  nanofluid 
concentration  moves  down  to  the  boundary  for 
increasing  values  of  <|). 

Fig.  14  displays  the  effects  of  thermal  radiation 
parameter  on  dimensionless  temperature  distribution. 
It  is  evident  that  the  temperature  of  both  nanofluids 
diminished  due  to  radiation  for  both  nanofluids.  The 
effect  of  temperature  ratio  parameter  on  temperature 
distribution  is  illustrated  in  fig.  15.  It  is  seen  that  the 
temperature  ratio  parameter  enrich  the  thermal 
boundary  layer  thickness  for  its  increasing  values. 

Various  values  of  Schmidt  number  on  concentration  is 
portrayed  in  fig.  16.  When  we  increase  the  values  of 


@  IJTSRD  |  Available  Online  @  www.ijtsrd.com  |  Volume  -  1  |  Issue  -  6  |  Sep  -  Oct  2017 


Page:  1306 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 


Sc,  the  nanoparticle  concentration  decreases  when 
solid  volume  fraction  <j>  =  0.05. 

Tables  3,  4  and  5  display  the  effect  of  various 
governing  physical  parameters  of  the  problem  on  the 
skin  friction  coefficient,  wall  concentration  and  > 
temperature  of  the  wall.  Table  3  presents  the  skin 
friction  coefficient  for  different  values  of  magnetic 
interaction  parameter,  suction  parameter,  shrinking 
parameter  and  solid  volume  fraction.  The  skin  friction 
coefficient  enhances  due  to  the  effect  of  magnetic  ^ 
field,  suction  and  solid  volume  fraction. 

The  numerical  values  of  wall  concentration  are  given 
in  table  4.  It  is  observed  that  the  shrinking  parameter 
enhances  the  concentration  of  the  wall.  In  the  mean 
while  magnetic  interaction  parameter,  solid  volume  > 
fraction,  suction  parameter  and  Schmidt  number 
diminishes  the  concentration  of  the  wall. 

From  table  5,  it  is  inferred  that  the  effect  of  magnetic 
field,  suction  and  thermal  radiation  over  the  wall 
temperature  is  to  suppress  it  while  shrinking 
parameter,  solid  volume  fraction  and  temperature  ^ 
ratio  parameter  enlarge  the  values  of  temperature  of 
the  wall. 

V.  CONCLUSIONS 

> 

The  numerical  solution  of  a  two  dimensional,  laminar, 
nonlinear,  steady  boundary  layer  flow  of  a  nanofluid  over^ 
an  incompressible  steady  shrinking  surface  subjected  to 
prescribed  heat  and  mass  flux  in  the  presence  of  thermal 
radiation  and  magnetic  field  with  suction  is  obtained  for 
various  values  of  governing  physical  parameters  in  this 
investigation. 


parameter,  solid  volume  fraction  and  temperature 
ratio  parameter  while  opposite  trend  is  occurred  in 
temperature  distribution  for  various  values  of 
magnetic  interaction  parameter,  suction  parameter 
and  thermal  radiation  parameter. 

An  increase  in  magnetic  interaction  parameter, 
suction  parameter,  solid  volume  fraction  and 
Schmidt  number  is  to  decrease  the  concentration 
distribution.  But  the  shrinking  parameter  raises  the 
nanofluids  concentration. 

The  important  physical  quantity  such  as  skin 
friction  coefficient  of  the  problem  increased  with 
increases  the  values  of  magnetic  interaction 
parameter,  suction  parameter  and  solid  volume 
fraction. 

The  temperature  at  the  wall  diminishes  when  we 
increase  the  physical  parameters  such  as  magnetic 
interaction  parameter,  suction  parameter  and 
thermal  radiation  parameter.  Whereas  the  wall 
temperature  boosts  for  increasing  the  shrinking 
parameter,  solid  volume  fraction  and  temperature 
ratio  parameter. 

An  increment  in  magnetic  field,  suction 
parameter,  solid  volume  fraction  and  Schmidt 
number  reduces  the  concentration  of  the  wall 
while  it  increases  with  rises  the  values  of 
shrinking  parameter. 

Ag  water  nanofluid  has  highest  skin  friction 
coefficient  as  compared  to  cu  water  nanofluid. 

The  numerical  values  of  wall  temperature  and 
concentration  of  the  wall  of  copper  water 
nanofluid  is  superior  to  silver  water  nanofluid. 


The  numerical  results  predict  the  effects  of  different 
physical  parameters  such  as  magnetic  interaction 
parameter,  shrinking  parameter,  suction  parameter, 
solid  volume  fraction,  radiation  parameter,  temperature 
ratio  parameter  and  Schmidt  number. 

The  main  results  of  the  work  can  be  summarized 
follows: 

>  The  momentum  boundary  layer  thickness 
enhanced  due  to  the  effect  of  magnetic  interaction 
parameter,  suction  parameter  and  solid  volume 
fraction.  Consequently  the  shrinking  parameter 
diminished  the  nanofluid  velocity  for  both  copper 
water  nanofluid  and  silver  water  nanofluid. 

>  The  temperature  distribution  in  the  boundary  layer 
flow  region  is  enriched  by  the  effect  of  shrinking 


TABLE  3:  VARIATIONS  IN  SKIN  FRICTION 
COEFFICIENT  FOR  VARIOUS  VALUES 
OF  M2,|3,(|),  AND  Fw  WHEN  N  =  3.0 


as 


Ejl 

P 

<|> 

f 

lw 

Cu  -  Water 

Ag  -  Water 

■ 

,  \l/2 

C  f  (Re,) 

/  \l/2 

C f  (Rex ) 

0.5 

1.5 

0.05 

2.5 

4.28818 

4.55550 

1.0 

4.57177 

4.82298 

1.5 

4.82269 

5.06254 

2.0 

5.05009 

5.28137 

2.0 

1.5 

0.0 

2.5 

3.90989 

3.90989 

0.05 

5.05009 

5.28137 

0.1 

6.24356 

6.73162 

0.2 

8.66328 

9.68743 

2.0 

1.5 

0.05 

2.0 

4.08432 

4.24979 

2.2 

4.46966 

4.66195 

2.3 

4.66285 

4.86827 

2.5 

5.05009 

5.28137 
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TABLE  4:  CONCENTRATION  OF  THE  WALL  ^(0)  FOR  VARIOUS  VALUES  OF  M2,p4,  Fw  AND  SC 


M2 

P 

<l> 

fw 

Sc 

Cu  -  Water 

Ag  -  Water 

<K°) 

m 

0.5 

1.5 

0.05 

2.5 

0.5 

0.89231 

0.88786 

1.0 

0.88770 

0.88390 

1.5 

0.88398 

0.88065 

2.0 

0.88088 

0.87789 

2.0 

1.5 

0.05 

2.5 

0.5 

0.88088 

0.87789 

2.0 

0.88128 

0.87910 

2.5 

0.88358 

0.88041 

3.0 

0.88512 

0.88184 

2.0 

1.5 

0.0 

2.5 

0.5 

0.89001 

0.89001 

0.05 

0.88088 

0.87789 

0.1 

0.87578 

0.87109 

0.2 

0.87365 

0.85836 

2.0 

1.5 

0.05 

2.0 

0.6 

1.16205 

1.15627 

2.2 

1.02963 

1.02524 

2.3 

0.97450 

0.97066 

2.5 

0.88088 

0.87789 

2.0 

1.5 

0.05 

2.5 

0.5 

0.88088 

0.87789 

0.6 

0.73106 

0.72879 

0.7 

0.62389 

0.62213 

0.8 

0.54364 

0.54224 

TABLE  5:  WALL  TEMPERATURE  0(0)  FOR  VARIOUS  VALUES  OF  M2,p,(|),ew,  NR  AND  Fw  WHEN  Pr  =  6.2 


M2 

P 

4> 

fw 

Nr 

0W 

Cu  -  Water 

Ag  -  Water 

0(0) 

0(0) 

0.5 

1.5 

0.05 

2.5 

2.0 

1.1 

0.15458 

0.15604 

1.0 

0.15383 

0.15536 

1.5 

0.15321 

0.15478 

2.0 

0.15268 

0.15428 

2.0 

1.5 

0.05 

2.5 

2.0 

1.1 

0.15268 

0.15428 

2.0 

0.16225 

0.16378 

2.5 

0.17420 

0.17557 

3.0 

0.18969 

0.19073 

2.0 

1.5 

0.0 

2.5 

2.0 

1.1 

0.12856 

0.12856 

0.05 

0.15268 

0.15428 

0.1 

0.18064 

0.18445 

0.2 

0.25325 

0.26442 

2.0 

1.5 

0.05 

2.0 

2.0 

1.1 

0.22240 

0.22423 

2.2 

0.18667 

0.18844 

2.3 

0.17347 

0.17518 

2.5 

0.15268 

0.15428 

2.0 

1.5 

0.05 

2.5 

1.0 

1.1 

0.23072 

0.23268 

2.0 

0.15268 

0.15428 

2.5 

0.13823 

0.13972 

3.0 

0.12882 

0.13023 

2.0 

1.5 

0.05 

2.5 

2.0 

0.8 

0.14422 

0.14566 

0.9 

0.14685 

0.14834 

1.1 

0.15268 

0.15428 

1.2 

0.15593 

0.15759 
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VI.  FIGURES 


Fig-^  Velocity  profdes  for  various  values  of  M2  Fig.3  Temperature  profile  for  various  values  of  M2 


Fig.4  concentration  profiles  for  various  values  of  M2  Fig.5  Velocity  profiles  for  various  values  of  (3 


Fig.6  Temperature  profiles  for  various  values  of  (3  Fig.7  Concentration  profiles  for  various  values  of  p 
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Fig.  8  Velocity  profiles  for  various  values  of  fw 


Fig.9  Temperature  profiles  for  various  values  of  fw 


Fig.  1 0  Concentration  profiles  for  various  values  of  fw  ^ig> 1 1  Velocity  profiles  for  various  values  of  (|) 


Fig.12  Temperature  profiles  for  various  values  of  4»  Fig.13  Concentration  profiles  for  various  values  of  § 


Fig.  1 4  Temperature  profiles  for  various  values  of  Nr  Fig.15  Temperature  profiles  for  various  values  of  0W 
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Fig.  16  Concentration  profiles  for  various  values  of  Sc 
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